Competition of ferromagnetism and superconductivity in Sc 3 InB 



B. WiendlochaQ J. Tobola, and S. Kaprzyk 
Faculty of Physics and Applied Computer Science, 
AGH University of Science and Technology, al. Mickiewicza 30, 30-059 Krakow, Poland 



in 
o 
o 



in 

^3 

o 

o 

Or 
c3 



^3 

o 

(N 
> 
m 

On 
O 

in 

o 

-i— > 
a 



i 

c 

o 
o 



D. Fruchart 

Laboratoire de Cristallographie, CNRS, BP 166, 



Grenoble Cedex 9, France 



J. Marcus 

Laboratoire d'Etude des Proprietes Electroniques des Solides, 
CNRS, BP 166, 38042 Grenoble Cedex 9, France 

We present results of electronic structure calculations for intermetallic perovskite ScalnB with 
Full-Potential KKR-LDA method. ScalnB is very promising candidate for a new superconductor 
(related to 8 K MgCNis) and can be regarded as a boron-inserted cubic Sc3ln, which is a high- 
pressure allotropic form of the hexagonal weak ferromagnet Scaln. We predict that cubic Scsln can 
be also magnetic, whereas Sc3lnB having large DOS in the vicinity of Ef exhibits non-magnetic 
ground state. Estimation of the electron-phonon coupling for Sc3inB gives A ~ 1. Furthermore, 
the effect of vacancy in Sc3lnBi_a; and antisite disorder in Sc3(In-B) on critical parameters is also 
discussed in view of KKR-CPA method. All theoretical results support possibility of the super- 
conductivity onset in ScsInB. Preliminary experimental measurements established the transition 
temperature close to 4.5 K, with a very abrupt change in susceptibility and a correlated drop of the 
resistivity when cooling down. 

PACS numbers: 74.10.+v,74.25. Jb,74.25.Kc,74.62.Dh,74.70.Ad,75.50.Ce 



I. INTRODUCTION 

Recent discovery of superconductivity in the inter- 
metallic perovskite MgCNi 3 [lj with T c ~ 8 K was a 
great surprise due to large Ni contents and this com- 
pound was rather expected to be near ferromagnetic crit- 
ical point 0- It was established that electron-phonon 
mechanism is responsible for superconductivity in this 
material, but few details are still not clear. NMR t\ 
relaxation time experiments 3'] or specific heat measure- 
ments [4| have resulted in typical behaviors, supporting 
s-wave type pairing with electron-phonon coupling con- 
stant A ~ 0.8. Conversely, rather unconventional be- 
haviors have been observed in other experiments (e.g. 
increase of critical temperature with pressure [jj or un- 
usual low-temperature behavior of London penetration 
depth Al(T) Q), which can be partly connected with 
spin fluctuations likely appearing due to vicinity of fer- 
romagnetism. Moreover, complex dynamical properties 
of this superconductor, as e.g. soft-mode behaviors and 
instability of Ni vibrations 0, make its theoretical 
analysis quite cumbersome. 

In this paper we report on predictions of superconduc- 
tivity in related ScsInB compound, which also seems to 
be close to magnetism limit due to weak ferromagnetic 
properties of SC3IJ1. It was already revealed that Sc3lnB 



'email: bartekwOfatcat. ftj.agh.edu. pi 

This work was presented at the European Conference Physics of 
Magnetism '05, 24-27.06.2005, Poznan, Poland. 



crystallizes in a perovskite structure (space group Pm- 
3m, CaTi0 3 type) with lattice constant a = 4.56 A [Io| . 
However, the same number of valence electrons of boron 
and indium elements may give rise to lattice instabilities 
and In/B antisite defects are plausible. 



II. THEORETICAL STUDY 

Electronic structure calculations were performed by 
Full Potential KKR method within the LDA 

framework employing von Barth-Hedin formula for the 
exchange-correlation potential. Figure ^ presents elec- 
tron density of states (DOS) in Sc 3 InB compound. The 
Fermi level (Ep) is located on the decreasing slope of 
large DOS peak with the n(Ep) value being as large 
as ~ 90 states/Ry. The main contributions can be at- 
tributed to Sc (d-states) and In (p-states) (see Tab. QJ. 
Noteworthy, appearance of a large DOS peak close to Ep, 
coming essentially from transition metal Ni-cf states, was 
also characteristic of electronic structure of MgCNi3 (Mg 
contribution was negligible) The SC3L1B case is 

not similar, since In plays more active role in formation 
of electronic states near Ep due to one electron more on 
p-shell. 

a. Predictions of superconductivity in ScalnB The 
electron-phonon coupling strength was estimated by 
calculating McMillan-Hopficld parameters [T3 | from 
Gaspari-Gyorfly formulas 0] - |16| within Rigid Muf- 
fin Tin Approximation (RMTA). The electron-phonon 
coupling constant was then calculated from relation: 
A = ^2i'Hi/' m i{ <jJ 'i)i where i corresponds to atoms in the 
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TABLE I: Properties of SC3I11B. n(E F ) is given in Ry 1 /atom, n in mRy/as 2 /atom, lo in meV, M(uj 2 ) in mRy/as 2 



Atom 


n(E F ) 


n s (E F ) 


n p (E F ) 


n d (E F ) 


n f (E F ) 


V 




Vpd 


Vdf 


V(" 2 ) 


M{uj 2 ) 


A 


Sc 


20.26 


0.06 


2.34 


16.98 


0.88 


18.6 


0.0 


4.7 


13.3 


18.5 


75.7 


0.74 


In 


18.28 


0.04 


16.72 


1.16 


0.32 


1.5 


0.0 


1.5 


0.0 


18.5 


193.4 


0.01 


B 


6.72 


0.02 


5.92 


0.42 


0.36 


34.3 


0.0 


34.3 


0.0 


58.3 


180.9 


0.19 
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FIG. 1: Total electron DOS in Sc 3 InB. 



FIG. 2: Phonon DOS in Sc 3 InB. 



unit cell. The (uif) parameter was derived from phonon 
DOS F (to) |22j |. computed for minimum-energy lattice 
constant a — 8.610 (1 ag = 0.529 A) within Density 
Functional Perturbation Theory, using the PWscf pack- 
age [13 ■ 

In order to take into account the markedly different 
masses of Sc, In and B atoms, phonon spectrum was an- 
alyzed while focusing on two parts: low frequency region 
with predominantly In and Sc modes (peak at 13 meV 
comes from flat acoustic In modes) and high frequency 
region with B vibrations (above 50 meV). Noteworthy, 
similar separation of phonon DOS was earlier observed 
experimentally in MgCNi3 Q- Values of (u)f) used next 
to estimate A, were calculated separately for both regions 
(with cut line at 50 meV): (uif) for Sc and In are taken 
to be equal, and represent the low-frequency part, while 
(lo 2 b ) was taken from high frequency part of the spectrum. 
Results and estimation of A^ are presented in Table [IJ to- 
tal electron-phonon coupling constant is A = 0.94. 

Using McMillan formula [p| for critical temperature 
and applying typical value of Coulomb pseudopotential 
Li* = 0.13, we got T c ~ 12 K (we used (w)/1.2 with 
(lo) = 19.5 meV instead of wrj/1.45 in McMillan formula 
in practical computations). Both A and T c values be- 
long to moderate regime of superconducting parameters 
within the RMTA model, and are higher than in MgCNi3. 

We can also notice advantageous trend of electron- 
phonon coupling in SC3IX1B if comparing with MgCNi3 
superconductor. Our KKR calculations for MgCNi3 
showed that Ni has the largest McMillan-Hopfield pa- 



rameter: r/jvi = 20 mRy/a| with respect to other atoms 
{rjc = 9 mRy/a|, rj Mg - negligible). The value of rj Sc 
in SC3L1B is similar to rjNi in MgCNi3 but rjB is over 
three times larger than rj c (both In and Mg contribu- 
tions are much smaller). In view of recent C isotope 
effect measurements [lg] (very large etc — 0.54 coeffi- 
cient) and bearing in mind a particular sensitivity of T c 
on the carbon concentration, one can conclude that both 
C and Ni sublattices are important in superconductivity 
of MgCNi3. In view of these criterion T c in SC3I11B can 
be expected higher than in MgCNi3. 

b. Magnetic properties of Sc^In It seems that in- 
stability towards magnetism may also be present in 
Sc3lnBi_ x compound. First, we have studied electronic 
structure of both allotropic phases of Scsln. The hexag- 
onal compound (Ni 3 Sn-type, a = 6.42 A, c = 5.18 A 
H3) is well-known weak itinerant ferromagnet, while 
the cubic compound (Cu3Au-type, a — 4.46 A, syn- 
thesized under high pressure [l(|) has not been yet in- 
vestigated to our knowledge. Present KKR calculations 
showed that both allotropic forms of Scsln should ex- 
hibit magnetic ground state supported by magnetic mo- 
ment on Sc atoms, i.e. 0.26 \ib (hexagonal phase) and 
0.27 fiB (cubic phase). We should remind that exper- 
imentally observed magnetic moment is much weaker 
(~ 0.05^b/Sc in hexagonal phase), as already underlined 
in the previous LAPW calculations 0]. Besides, the 
superconductor-to-ferromagnet transition can appear in 
ScsInBi-^ if varying B content (KKR-CPA computa- 
tions are in progress). Thus, as already suggested in 
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MgCNi3 0] ! one can expect that the proximity of fer- 
romagnetic quantum critical point (resulting in enhanced 
spin fluctuations) may probably compete with supercon- 
ductivity. 



III. EXPERIMENTAL ANALYSES 

A ~1.5 g sample was first prepared by arc melting 
under high purity argon atmosphere (99.9995) start- 
ing from the appropriate proportions of the elements 
(purity > 99.95) to obtain the Sc3inB formula. The 
resulting small ingot was melted several times in or- 
der to insure homogeneity. Pieces of the ingots were 
made by using a metal mortar, and inspection by us- 
ing an optical microscope reveals the bright and homoge- 
neous aspect of the fractured surfaces. Then, XRD pat- 
terns were recorded at \x a (Cu) using a Bragg-Brcntano 
diffractometer equipped with a backscattering pyrolitic 
graphite monochromator. The diffraction pattern reveals 
the presence of a dominant amount of cubic phase with 
the addition of a minor impurity. Probably because no 
annealing procedure was applied, the crystallised state 
of the sample was not of the best and no effective crys- 
tal structure refinement was applied. However using the 
PowderCell code the two main phases were clearly 
identified and a rough determination of the cell param- 
eters was made possible. The main phase (~70 % vol.) 
is simple cubic of perovskite type of structure and the 
second one (~30 % vol.) is hexagonal of Ni2ln type of 
structure (space group PQ^/mmc). Accounting for the 
large difference in between the scattering lengths of the 
p-elements B and In, an estimate of the composition of 
the main phase was made. All the results are displayed 
on Table El 



TABLE II: Structure analysis of the material synthesized with 
the nominal composition ScalnB. 



Compound 


Type 


Cell parameters [A] 


Composition 


ScalnB 








(~ 70% vol.) 


CaTiOa 


4.66 4.66 4.66 


Sc 3 Ini. 3 B .7 


Sc2ln 








(~ 30% vol.) 


Ni 2 In 


5.05 5.05 6.30 


Sc 2 In 



The best agreement for the main phase composition 
leads to consider that the p-elements B and In are nei- 
ther fully ordered nor randomly distributed, with about 
0.70 In and 0.30 B atom on the la site, then 0.65 In and 
0.35 B atom on the lb site of the space group Pm-3m. 
Consideration for such a type of atomic disordering was 
already reported in literature [2l]]. Both the composi- 
tions of the two compounds (3-1-1 and 2-1) confirm that 
boron is a difficult element to combine, thus the remain- 
ing boron should be detected in the XRD pattern, unless 
it is difficult to evidence effectively as a very light and 
often poorly crystallised element. 



Susceptibility and resistivity measurements were made 
using a a.c. susceptometer in temperature ranging from 
50 to 2 K. Millimeter sized pieces of the ingot were mea- 
sured successively and lead to the same results. A typical 
record is shown on Figure^ thus revealing the onset of a 
superconducting state down to 4.4 K, the transition being 
very sharp with no detectable hysteresis loop. A change 
in the resistivity trace was also observed simultaneously. 



FIG. 3: A.c. susceptibility trace revealing the transition to 
superconducting state. 
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Unexpectedly, the rest of the ingot and the small pieces 
rapidly change their brittle and bright aspects after sev- 
eral hours let in ambient atmosphere. Besides, a tenta- 
tive to melt twice the main parts of the initial ingot was 
made, but at this time no evidence for any supercon- 
ducting transition was found. Then, several new synthe- 
ses were undertaken, thus operating as possible similarly 
as for the first one procedure. Again, the new ingots 
we obtained do not display any transition down to 2 K. 
As revealed by the theoretical derivations (see below), 
the superconducting state and the related transition look 
fairly dependant on the p-element ordering and the sto- 
ichiometry in the la and lb sites. So, we anticipate that 
during and after the different melts, the new samples do 
not exhibit the same atom ordering as it was resulting 
from the first attempt. New syntheses are now scheduled 
to be undertaken using different methods and techniques 
in order to achieve optimized compounds. 

Because of experimental problems with synthesis, ef- 
fect of boron site vacancy on electronic structure near 
Ep and superconducting properties was simulated us- 
ing KKR method with coherent potential approximation 
(CPA). We have found significant change in r\ values 
not only for B, but also for Sc, i.e. for 7% boron de- 
ficiency (Sc3inBo.g3), employing the same values of (tof), 
total A decreased over 30% to 0.62, which resulted in 
T c ~ 4 K. In ScsInBo.ss coupling constant is so small 
(A ~ 0.4), T c < 0.5 K seems to be below the standard 
low-temperature measurements. Similar decrease (but 
less rapid than in the case of B vacancy) of critical pa- 
rameters was detected from KKR-CPA analysis when an- 
tisite In/B disorder increased. 
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